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Abstract 

Introduction: Circulating tumor cells (CTCs) that present mesenchymal phenotypes can escape standard methods of 
isolation, thus limiting possibilities for their characterization. Whereas mesenchymal CTCs are considered to be more 
malignant than epithelial CTCs, factors responsible for this aggressiveness have not been thoroughly defined. This study 
analyzed the molecular profile related to metastasis formation potential of CTC-enriched blood fractions obtained by 
marker unbiased isolation from breast cancer patients without (N— ) and with lymph nodes metastases (N+). 

Materials and Methods: Blood samples drawn from 117 patients with early-stage breast cancer were enriched for CTCs 
using density gradient centrifugation and negative selection with anti-CD45 covered magnetic particles. In the resulting 
CTC-enriched blood fractions, expression of CK19, MGBl, VIM, TWISTI, SNAIL, SLUG, HER2, CXCR4 and uPAR was analyzed with 
qPCR. Results were correlated with patients' clinicopathological data. 

Results:CTCs (defined as expression of either CK19, MGBl or HER2) were detected in 41% (20/49) of N- and 69% (34/49) of 
N+ patients (P = 0.004). CTC-enriched blood fractions of N+ patients were more frequently VIM (P = 0.02), SNAIL (P = 0.059) 
and uP/4/?-positive (P=0.03). Positive VIM, CXCR4 and uPAR status correlated with >3 lymph nodes involved (P = 0.003, 
P = 0.01 and P = 0.045, respectively). In the multivariate logistic regression MGBl and WM-positivity were independently 
related to lymph node involvement with corresponding overall risk of 3.2 and 4.2. IVloreover, mesenchymal CTC-enriched 
blood fractions (CK19-A/IM+ and MGB1+ or HER2+) had 4.88 and 7.85-times elevated expression of CXCR4 and uPAR, 
respectively, compared with epithelial CTC-enriched blood fractions (CK19+A/IM- and MGB1+ or HER2+). 

Conclusions: Tumors of N-i- patients have superior CTC-seeding and metastatic potential compared with N- patients. These 
differences can be attributed to VIM, uPAR and CXCR4 expression, which endow tumor cells with particularly malignant 
phenotypes. 
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Introduction 

Presence of circulating tumor cells (CTCs) in blood of patients 
with epithelial cancer was first demonstrated by Ashworth in 1 869 
[1]. As techniques were developed to capture, enumerate and 
characterize circulating and disseminated tumor cells, significant 
progress was made in understanding metastatic processes [2,3,4,5]. 
The number of CTCs detected in blood samples carry prognostic 
information in early [6,7] and metastatic breast cancer [8,9]. Also, 
CTCs detected via PCR-based methods (without the possibility of 
cell enumeration) have been associated with poor prognosis in a 
number of studies [10,11], described in a recent meta-analysis 
[12]. 



As CTCs originate from the epithelium, use of epithelial 
markers (eg, cytokeratins, EpCAM) for their detection seems 
reasonable. Cytokeratin 19 (CK19) is a cytoskeletal protein of 
epithelial cells (both normal and cancerous) and is widely used for 
detection of CTCs [13,14,15] and DTCs [15,16,17,18]. However, 
discovery of epithelial-mesenchymal transition (EMT) in cancer 
educes a reconsideration of CTCs as having exclusively epithelial 
phenotype [2,19,20,21]. Moreover, the role of EMT in cancer 
implies that detection methods that rely solely on epithelial 
markers (or other markers downregulated during EMT) are likely 
to miss the most aggressive fraction of CTCs [2,19,22]. Thus, to 
increase sensitivity it is suggested to include additional mammary 
transcripts, like mammaglobin 1 (MGBl), which was shown to be 
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a useful marker for detecting disseminated breast cancer cells iu 
blood [23,24,25], bone marrow [26,27] and lymph nodes 
[28,29,30]. Additionally to CK19 and MGBl, detection of HER2 
transcripts, which is frequently overexpressed in breast cancers, 
strengthened prognostic value of the RT-qPCR based CTCs 
detecting assay [23]. 

Activation of EMT is linked to motility, stem cell characteristics, 
enhanced chemo- and radiotherapy resistance [19,31,32,33]. The 
fraction of CTCs with a mesenchymal phenotype reportedly 
reaches almost 100% in the blood of some breast cancer patients 
[3]. Moreover, in some patients, disease progression during 
treatment was related to increased number of mesenchymal CTCs 
compared with their pre-treatment state [2] . The ability of tumor 
cells to metastasize can be modified by expression of various 
invasion and metastasis-related factors. Plasminogen activator, 
urokinase receptor (uPAR) constituting a part of uPA-PAI 
extracellular matrix degradation system might facilitating tumor 
cells invasion, migration and growth [34,35], uPAR was also shown 
to be amplified together with HER2 in breast cancer CTCs [36] 
and decreased expression of uPAR related to tumor cell dormancy 
[35]. Yet another protein, CXCR4 chemokine receptor, apart 
from being involved in metastases formation and migration of 
cancer cells to specific organs [37,38] is functionally linked with 
HER2 signalling and malignant progression. CXCR4 expression 
is enhanced by HER2, which can together act in multiple steps of 
metastatic cascade [39]. 

Inherent increased malignancy of mesenchymal CTCs could 
also contribute to higher metastatic potential, which in early-stage 
breast cancer could be measured by lymph-node involvement. We 
have hypothesized that CTCs isolated from lymph node-negative 
(N— ) and -positive (N+) patients could differ in expression of 
malignancy-associated genes. We therefore used a marker- 
unbiased CTC-enrichment method that enriches both epithelial 
and mesenchymal CTCs, in which we measured expression of 
mammary epitheUal transcripts (CK19, MGBl), EMT-related 
factors {VIM, TWISTl, SNAIL, SLUG) and invasion- and 
metastasis- related genes {HER2, CXCR4, uPARj. 

Methods 

Patients 

The study included 1 1 7 breast cancer patients, stages I-III, 
treated in the Medical University Hospital in Gdansk between 
April 201 1 and May 2013. Tumor stage and node positivity were 
defined according to AJCC cancer staging manual version 7. 
Micrometastases were considered and there was one case in 
examined group identified, which we classified as N+ patient. 
There were no patients with isolated tumor cells in the N- group. 
Cancers were graded according to modified Bloom and Richard- 
son system based on semi-quantitative method for assessing 
histological grade in breast tumours [40]. Median age of the 
patients \vas 61 years (28-89 years) (Table 1). Inclusion criteria 
were primary operable breast cancer confirmed by histological 
examination, and signed consent form. Peripheral blood samples 
(5-10 mL) were drawn to EDTA-coated tubes before tumor 
excision and therapy initiation. The first few milliliters of blood 
were discarded to minimize possibility of keratinocyte contamina- 
tion. Samples were stored at 4 °C until analysis, but no longer than 
24 h. Blood samples from 1 7 healthy women (median age 36 
years; range 20-73 years) and three women with benign breast 
disease (ductal carcinoma in situ - DCIS, median age 65 years; 
range 46-7 1 ) were similarly drawn and processed. 

Median follow-up time was 1.5 years (0.2 to 2.2 years). 
Including the last foUow-up data six deaths were observed, which 



is insufficient for performing survival analysis; however, follow-up 
data continue to be collected. 

Hormone receptors status (ER and PgR) was evaluated using 
classical Allred scoring method with 3 being a cut-ofF point for 
positive result. Standard criteria for evaluating HER2 positivity 
were applied, being 3-1- score in immunohistochemistry or positive 
result in fluorescence in situ hybridization (FISH), as previously 
described [41]. 

Ethics Statement 

The study was approved by the local Ethical Committee of the 
Medical University of Gdansk and the manuscript was prepared 
according to the REMARK criteria [42]. 

CTCs isolation/enrichment 

For the immunofluorescence experiment, blood samples were 
first centrifuged at 200 g for 20 min at 20 °C to remove excessive 
platelets that hampered cell adhesion to polylysrne slides after 
CTC-enrichment. The top serum layer, which contained platelets, 
was collected and discarded; the remaining fraction was then 
processed with the full blood sample intended for RNA isolation. 
Blood samples, or platelet-removed fractions (5 mL), were 
subjected to CTC-enrichment as described before [43] 
(Figure 1). Briefly, phosphate buffered saline (PBS) diluted blood 
sample was transferred into a 15 ml tube containing two-layer 
density gradient (upper and lower gradient) and centrifuged to 
separate tumor cells containing fraction from erythrocytes and 
blood cells. The tumor cell containing fraction was collected and 
subjected to further depletion of CD45-positive cells with anti- 
CD45-covered magnetic particles (CD45 Dynabeads, Invitrogen). 
After depletion, the obtained cells pellets were used in the RNA 
isolation procedure, or if cells were to be visualized by 
immunofluorescence, pellets were suspended in 1 mL PBS buffer 
and spun down on 2-4 polylysine-coated glass slides in Rotofrx 
32A (Hettich). 

To evaluate cell recovery rates in the developed CTC- 
enrichment method, a spike-in experiment was carried out twice. 
On average, 6 and 11 MDA-MB-361 cells (5 and 11 in the first 
experiment; 6 and 10 cells in the second experiment) were spiked 
to peripheral blood samples from healthy women, which were then 
processed as the patients' samples. Four cytospins were prepared 
from each sample, slides were stained for CK19 and CD45 as in 
the immunofluorescence experiment and analyzed under Axiovert 
200 (Zeiss) microscope. 

RNA isolation and qRT-PCR 

RNA was isolated with Trizol reagent (Invitrogen) according to 
manufacturer's instructions. Up to 1 0 nl of RNA was used in a 
reverse transcription reaction (Transcriptor FirstStrand Synthesis 
Kit, Roche) with random hexamers according to manufacturer's 
instructions. Presence of inhibitors in the isolated RNA was tested 
using exogenous RNA molecules (Solaris RNA Spike Control Kit, 
Thermo Scientific) according to manufacturer's instructions, 
adding 1 |.il of lOOx Solaris Spike solution per 1 (xg of RNA. 
cDNA was diluted and 10 ng of cDNA (4 (tL) was used in a single 
20-pL qPCR reaction. Universal PGR Mastermix (with UNG 
AmpErase) and TaqMan Gene Expression Assays were used to 
measure expression of nine genes of interest: TWISTl 
(Hs00361186_ml; UniGene Hs.644998), ShAIL (also known as 
SJVAIl, Hs00195591_ml; UniGene Hs.48029), SLUG {aho known 
as Smi2, Hs00950344_ml; UniGene Hs.360174), CK19 
(Hs01051611_gH; UniGene Hs.654568), MGBl 

(Hs00935948_ml; UniGene Hs.46452), HER2 

(Hs99999005_mH; UniGene Hs.446352), CXCR4 
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Table 1. Patients' characteristics (N = 117). 



Variable 


Number of cases (%) 


Age - median (range) 


61 


(28-89) 


T stage 


Tl 


50 


(43) 


T2 


58 


(50) 


T3 


5 


(4) 


T4 


3 


(3) 


Missing data 


1 


(1) 


N stage 


N- 


60 


(51) 


N+ 


57 


(49) 


Grade 


Gl 


15 


(13) 


G2 


63 


(54) 


G3 


39 


(33) 


HER2 status 


Negative 


88 


(75) 


Positive 


26 


(22) 


Missing data 


3 


(3) 


ER status 


Negative 


23 


(20) 


Positive 


94 


(80) 


PR status 


Negative 


29 


(25) 


Positive 


88 


(75) 


Histological type 


Ductal 


88 


(75) 


Lobular 


16 


(14) 


Other 


12 


(10) 


Missing data 


1 


(1) 


doi:l 0.1 371 /journal.pone.0093901 .tOOl 



(Hs00237052_ml; UniGene Hs.593413), uPAR (Hs00182181_ml; 
UniGene Hs.466871), VIM (Hs00185584_ml; UniGene 
Hs.455493) and two reference genes GAPDH (Hs99999905_ml; 
UniGene Hs.544577) and YWHAZ (Hs03044281_gl; UniGene 
Hs. 492407). Reference genes were chosen for their expression 
stability as assessed by geNorm according to the provided manual 
[43]. qPCR reactions were performed in duplicate on 96-well 
plates in a CFX96 thermal cycler (Bio-Rad). Cycling parameters 
were as follows: 10 min 95 °C, 45 cycles of 1 min at 60 °C 
followed by 30 s at 95 °C. For the TaqMan assays detecting 
genomic DNA [CK19 and TWHA^, sample-specific controls 
(10 ng of untranscribed RNA) were included on every plate. 
Samples were assumed to be genomic DNA-free if diflFerences 
between control and test samples were a5 CT. Gene expression 
was calculated using a modified AACt method that corrects for a 
run-to-run variation [44]. Therefore, every plate included an 
inter-run calibrator that allowed calculation of calibrator-normal- 
ized relative quantities using qBasePLUS software version 2.1 
[44] . Gene expression levels were scaled to the minimal expression 
level of each gene in the patients' samples. 



Immunofluorescence 

As PCR-based techniques do not allow for cells visualization, 
immunofluorescence experiments were carried out for some 
patient samples to visuahze CTCs isolated using the developed 
method (described in "CTCs isolation/enrichment"). Double- 
staining experiments were performed for the following proteins 
combinations: CK19 and CD45, CK19 and MGBl, CXCR4 and 
CD45, HER2 and CD45, SNAIL and CD45. Primaiy mouse 
monoclonal (CK19, CXCR4, SNAIL, HER2) or rabbit polyclonal 
(CD45, MGBl) antibodies were used (Santa Cruz Biotechnology; 
Abnova for anti-SNAIL antibody). Cytospins, prepared on 
polylysine-coated shdes as described above, were fixed with cold 
methanol for 5 minutes and incubated with appropriate pairs of 
primary antibodies for 1 hour at room temperature. Pairs of 
mouse monoclonal antibodies and rabbit polyclonal antibodies 
were used. Slides were then stained for 1 h with corresponding 
secondary antibodies: Sheep Anti-Mouse DyLight 549 and Goat 
Anti-Rabbit DyLight 488 JIR). Cell nuclei in all slides were 
stained with 4,6-diamidino-2-phenylindole (Sigma). Slides were 
evaluated using inverted fluorescent microscope Axiovert 200 
(Zeiss) and analyzed with Axio Vision software (Zeiss). 

Statistical analysis 

Statistical analyses were performed with STATISTICA soft- 
ware version 10. Categorical variables were analysed using 
contingency tables with x2 statistics or Fisher's exact test where 
applicable. Spearman rank coefficient was used when continuous 
variables were correlated; Mann-Whitney test for analyzing 
continuous variables distribution in two groups. Logistic regression 
analysis was used to identify the gene predictors of lymph node 
involvement. Univariate predictors significant with a value of p^ 
0.10 were entered into a step-wise backwards multivariate model 
to identify those with independent prognostic information. 

Unsupervised hierarchical clustering analysis was performed 
using GenePattern version 6 software (http://genepattern. 
broadinstitute.org/) [45]. Spearman's rank correlation was applied 
for columns and rows distance measure. Pairwise complete-linkage 
was used as a clustering method. Samples with missing gene 
expression values were removed from the analysis. 

Results 

CTC-enriched blood samples characterization 

For the average 6 and 11 MDA-MB-361 cells added to blood 
samples, the CTC-enrichment method showed average recovery 
rate of 54% and 72%, respectively. In every case, 2-4 cells were 
lost during sample processing. Double staining (staining of CD45 
and another protein of interest) revealed CK19+/CD45 — , 
HER2+/CD45-, CXCR4+/CD45- and SNAIL+/CD45- cells 
isolated from breast cancer patients (Figure SI). In CK19/MGB1 
double staining, following phenotypes occurred: CK19+/MGB1+, 
CK19+/MGB1- and CK19-/MGB1+ (Figure SI). 

Functional RNA was successively isolated from 84% (98/ 1 1 7) of 
patient samples and 65% (11/17) of controls. None of the control 
samples was positive for TWISTl, SMIL, SLUG, HER2 or uPAR; I 
out of 1 1 showed CXCR4 expression (relative gene expression level 
4.54); and 9 (82%; 9/11) showed EH^f expression (median relative 
gene expression level 11.44, range 0-30.67). Expression oi MGBl 
and CK19 was analyzed in a qualitative manner. Patient samples 
were considered positive when measured relative gene expression 
level was higher than the highest measured relative expression 
level in the control samples. Table 2 presents the number of 
samples expressing analyzed genes according to this positivity 
criterion. Fifty-five percent (54/98) of the patient samples with 
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Blood sample Before After 

collection Density gradient centrifugation 



Negative selection 
with anti-CD45 
magnetic beads 



RNA isolation & 
reverse transcription 



Gene expression analysis 
with qPCR 



Figure 1. Schematic representation of blood samples analysis process. Collected blood sample was layered on density gradient. After 
centrifugation fraction containing tumor cells (marked by green box) was collected and subjected to negative selection. After CD45-depletetlon, CTC- 
enrlched blood sample was subjected to RNA Isolation, reverse transcription and gene expression analysis with qPCR. 
dol:10.1371/journal.pone.0093901.g001 



invasive tumors were positive for either CK19, MGBl or HER2. In 
total, 26% of the samples expressed at least one EMT marker 
{VIM, TWISTl, SNAIL or SLUG). No healthy control expressed 
CTCs markers, whereas in all three cases of DCIS expre.ssion of at 
least one CTC marker was detected (Table 2 and Figure S2). 

CTC-enriched samples expressing at least one of CTCs markers 
— MGBl or HER2 were divided by CK19 and VIM gene expression 
status into four group.s with positive or negative expression status. 
Phenotype frequencies were as follow: CK19+/VIM— 17% (16/ 
94), CK19+/VIM+ 1% (1/94), CK19-/VIM+ 13% (12/94) and 
CK19-/VIM- 69% (65/94). CK19+/VIM+ sample was excluded 
from statistical analysis to avoid biased classification to either 
CK19+/VIM- or CK19-/VIM+ group. Expression oi CXCR4 and 
uPAR was higher in F/M-positive fraction {CK19- /VIM+) than in 
the CK19-/VIM- CTC-enriched blood fraction (Figure 2). In 
CK19—/VIM+ CTC-enriched blood fractions, average expression 
of CXCR4 and uPAR was 4.88 and 7.85-fold higher than in 
CK19-/VIM- fraction, respectively {P < 0.00001 for both). 

Correlations between gene expression of examined 
markers 

SNAIL expression was found in 8% (8/96) of the samples; its 
expression correlated with positive status of VIM {P= 0.03), HER2 
{P= 0.0002), CXCR4 {P= 0.006) and uPAR (P= 0.02) (Table 3). 
CTC-enriched //£ft2-positive blood samples were more frequent- 
ly F/M-positive (38% vs. 12%, P= 0.004), &¥Att-positive (24% vs. 
0%, P= 0.0002), OrCft^-positive (53% vs. 33%, P=0.05) and 
wP^-positive (100% vs. 34%, P<0.00001) than HER2-negative 
samples (Table 3). Moreover, 82% (32/39) of CTC-enriched 
CXCS^-positive samples were also wiMP -positive (P= 0.0001; 
Table 3). Spearman rank correlations between relative expression 
of the tested genes showing similar gene dependence are presented 
in Table SI. 



Correlation with clinicopathological data 

In CTC-enriched blood fractions of N+ patients, at least one of 
the markers - CK19, MGBl or IIER2 - was more frequently (69%, 
34/49) detected than in N- patients (41%, 20/49; P= 0.004, 
Table 4). Concerning phenotype of CTC-enriched blood 
fractions, there was no difference between lymph node involve- 
ment in patients with epithelial and mesenchymal CTC-enriched 
blood fraction phenotype (P— 0.69), however mesenchymal 
phenotype was more frequendy found in patients with more than 
3 lymph nodes involved —50% in comparison to patients with 
epithelial CTC-enriched blood fraction phenotype —7% 
(P= 0.008, Table 5). Moreover, in N+ patients, CTC-enriched 
blood fractions were more often MGBl, SNAIL, VIM and uPAR- 
positive {P= 0.03, P= 0.059, P= 0.02 and P= 0.03, respectively) 
(Table 5). Positive status of CXCR4 was associated with higher T 
stage [P= 0.01). There was no correlation between HER2 status of 
primary tumors (PT) and IIER2-stMm of CTC-enriched blood 
fractions (P = 0.87) (Table 5). 

Multivariate analysis revealed that independent predictors of 
lymph nodes involvement were positive status oi MGBl - OR 3.2 
(95% CI 1.1-9.2, P= 0.029) and FiM- OR 4.2 (95% CI 1.3-13.5, 
P=0.01) (Table 4). 

In hierarchical clustering, the study population was divided into 
two main groups, which differed in expression of VIAL, CXCR4, 
uPAR, HER2 (Figure 3). Patients in the cluster with elevated 
expression of these genes showed more frequent lymph node 
involvement (58%) than patients from the cluster with lower 
expression (35%; P=0.03). 

Discussion 

Dissemination of cancer cells is an early event and disseminated 
tumor cells can be found in bone marrow of patients with 
carcinoma in situ [46]. Nevertheless, in metastatic breast cancer 
more CTCs are seen than in early breast cancer patients 
[13,16,47]. No clear association is apparent between lymph node 
involvement and CTC detection rate. Some studies show similar 
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Figure 2. Relative expression ievei of CXCR4 and uPAR in CTC-enriclied biood fractions. CTC-enriched blood fractions positive for either 
MGB1 or HER2 were divided according to VIM and CK19 expression status Into three groups {CK19-/VIM—, CK19+A/IM—, CK19-/VIM+). Error bars 
depict standard error. * - statistically significant difference In CXCR4 and uPAR relative expression level (P<0.00001 for both) between CK19—A/IM— 
and CK19-/VIM+; ** - statistically significant difference In CXCR4 (P=0.00002) and uPAR (P = 0.00004) relative expression level between CK19^A/IM- 
and CK19-/VIM+. 
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CTCs detection rate in N— and N+ breast cancer patients [23,48], 
which would indicate similar tumor seeding potential, but 
dissimilar colonization potential for disseminated cells. However, 
preliminary results from a large SUCCESS trial revealed a 
correlation between CTC presence and lymph node involvement 
[49], implying different seeding, and possibly colonization, 
potentials in N— and N+ patients. To explore differences in 
seeding and colonization potential we analyzed expression of 
mammary epithelial transcripts [CK19, MGBl), EMT-related 
factors [TWISTl, SNAIL, SLUG, VL\T) and invasion and metasta- 
sis-related genes (HER2, CXCR4, uPAR} in CTC-enriched blood 
fractions from N- and N+ breast cancer patients. Considering the 
unresolved discussion concerning markers for CTC detection and 
characterization [50] and drawbacks of CTC-positive selection, we 
decided to apply a negative selection method for CTC- 
enrichment. Similar to prior research identifying CTCs with 
PCR-based methods, we used a "classical" definition of CTCs 
presence, meaning expression of either CK19, MGBl or HER2 
[23]. However, growing body of evidence shows that mesenchy- 
mal CTCs are not a rare event in cancers, what made us apply an 
extended definition of CTCs phenot^pes - epithelial or mesen- 
chymal - which includes VIM as an additional marker related to 
mesenchymal state. As a result we consider CTCs-enriched blood 
fractions being CK19+/VIM- and MGB1+ or HER2+ as 
containing CTCs in epithelial state, and CK19— /VIM+ and 
MGB1+ or IIER2+ fractions as carrying mesenchymal CTCs. 
Even though the phrase describing MGBl as "mammary epithelial 
transcript" implies it is present in mammary epithelium, it does 
not mean that it is expressed only by cells in the epithelial state. 
Although MGBl function is not known, structurally related 
proteins are a group of secretory proteins binding steroid ligands 
that might present anti inflammatory activity [51,52]; thus 



implying related function in human breast tissue, not restricted 
to epithelial state. 

In our study, we found expression of either OKI 9, MGBl or 
HER2 in 55% of the patients, which agrees with previous studies of 
early breast cancer patients [53,54,55]. Unlike in healthy controls, 
in DCIS cases expression of at least one CTC marker was 
detected, what supports results of the early cancer dissemination 
model [46]. We have observed significanfly different CTCs 
detection rate (defined as expression of either OKI 9, MGBl or 
HER2) in N- and N+ patients (41% vs. 69%, respectively, 
P= 0.005; Table 4), which supports the notion that tumors from 
N+ patients have superior seeding potential compared with N— 
patients. We would not have drawn the same conclusion if CK19 
was the only CTC marker in our study, as in that case there would 
be no significant difference in CTC marker detection rate between 
N- and N+ patients (20% vs. 33%,, P=0.17; Table 4). Also, the 
study of Pecot indicates that cytokeratin is not an ultimate marker 
for CTC identification as cytokeratin-negative cancer cells can be 
found both in circulation and within FTs [56]. Multimarker 
approach for CTCs detection is therefore a necessity [14,57]. Our 
immunofluorescence results also showed CK19— /MGB1+ cells in 
CTC-enriched blood fractions, which was also seen in prior 
research [58]. Analysis of expression of EMT markers, and of 
invasion- and metastasis-related genes revealed that CTC- 
enriched blood fractions of N-l- patients are more frequentiy 
FZ'Vf-positive and a/MiJ-positive (and show a trend toward being 
SMAIL-positive); and patients with more than three involved lymph 
nodes are often VIM, CXCR4 and wiMiJ-positive. Increased 
metastatic abilities of tumor cells in N+ patients could be therefore 
attributed to expression of CXCR4, uPAR and VIM (and possibly 
SNAIL). CXCR4 is a receptor for CXCL12 chemokine, which is 
secreted by the common sites of breast cancer metastasis, including 
lymph nodes [3]. Thus, CXCR4 expression on CTCs could 
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Table 4. Results of univariate and multivariate logistic regression analysis in relation to lymph node involvement. 





Gene expression status 


Number of cases (%) 


Univariate analysis 




Multivariate analysis 






N- 


N+ 


OR (95% CI) 


P 


OR (95% CI) 


P 


CKig 


Negative 


39 (80) 


33 (67) 


1 .9 (0.75-4.8) 


0.17 






Positive 


10 (20) 


16 (33) 










MGB1 


Negative 


40 (85) 


31 (66) 


2.9 (1.06-8.2) 


0.03 


3.2 (1.1-9.2) 


0.029 


Positive 


7 (15) 


16 (34) 










HER2 


Negative 


32 (67) 


29 (62) 


1 .24 (0.5-2.9) 


0.6 






Positive 


16 (33) 


18 (38) 










VIM 


Negative 


43 (90) 


32 (70) 


3.76 (1.2-11.7) 


0.02 


4.2 (1.3-13.5) 


0.01 


Positive 


5(10) 


14 (30) 










TWISTl 


Negative 


47 (98) 


47 (98) 


1 (0.6-17) 


1 






Positive 


1 (2) 


1 (2) 










SNAIL 


Negative 


47 (98) 


41 (85) 


8.02 (0.92-69.9) 


0.06 




NS 


Positive 


1 (2) 


7(15) 










CXCR4 


Negative 


32 (65) 


25 (52) 


1.73 (0.76-3.9) 


0.19 






Positive 


17 (35) 


23 (48) 










uPAR 


Negative 


25 (52) 


15 (31) 


2.46 (1.06-5.7) 


0.03 




NS 


Positive 


23 (48) 


34 (69) 










CK19/MGB1/HER2 


Negative 


29 (59) 


15 (31) 


3.29 (1.41-7.64) 


0.005 


Not applicable 




Positive 


20 (41) 


34 (69) 











N- lymph node involvement absent, N+ lymph node involvement present. OR-overall risk; CI - confidence interval, NS-not statistically significant. 
doi:1 0.1 371 /journal.pone.0093901 .t004 



mediate homing of tumor cells to lymph nodes. Also, the role of 
uPAR in extracellular matrix degradation complex, migration and 
invasion could explain its metastasis-promoting function [34]. 
Amplification and increased expression of uPAR was observed 
before in CTCs [36], PTs, and disseminating tumor cells in bone 
marrow and in lymph nodes of breast cancer patients [59] . Co- 
expression of uPAR and HER2 shown by Meng et al. [36] also 
supports our observations in /ffiK2-positive CTC-enriched blood 
fractions, which were all M/L-Lft-positive. 

Although induction and maintenance of EMT are thought to 
require TlVISTl, SMAIL and SLUG, their expression was not 
found in CTCs expressing mesenchymal markers [2,60], nor were 
they highly expressed in our CTC-enriched blood fractions; only 
8% and 2% of samples were positive for either SMAIL or TWISTl, 
respectively and SLUG was not detected in any sample. Similar 
rates were detected by Mego et al., who analyzed CTC-enriched 
blood fractions depleted of EpCAM-positive CTCs [22]. Despite 
the small number of >SAC4/Z-positive samples, SMAIL expression 
correlated with positive status of VIM, CXCR4, uPAR and IIER2. 
We also saw positive VIM status in 20% of the CTC-enriched 
blood fractions, which correlated with SMAIL, CXCR4 and HER2 
expression. This result supports data showing increased malig- 



nancy of CTCs with mesenchymal phenotype, possibly resulting 
from EMT [20,33,61]. No significant increase in expression of 
CXCR4 and uPAR in OKI 9+ /VIM- CTC-enriched blood fraction 
would speak for the passive model of dissemination, in which 
tumor cells are physically translocated into the vasculature (or the 
neovasculature formed around tumor cells) [62]. In that case, cells 
would not need a migratory phenotype, typically associated with 
EMT, and could remain in their epithelial phenotype, which 
according to several theories would strengthen their metastatic 
ability because of increased adhesive properties. We presume that 
CTCs with mesenchymal phenotype (lacking epithelial adhesion 
molecules) acquire additional malignant properties (expression of 
CXCR4 and uPAR) that allow successful lymph node colonization. 
It would be interesting to check whether expression of CXCR4 and 
uPAR are related to any specific metastatic pattern other than the 
lymph node metastasis seen in our study. 

Limitations of our study include relatively small sample size and 
short follow-up period, which hamper survival analysis. Moreover, 
the associations between variables are deduced from their inter- 
correlation, which does not necessarily inform about causal 
relationship. Therefore, our results should be seen as hypothesis- 
generating discoveries, and would obtain additional strength when 
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Figure 3. Unsupervised hierarchical clustering of relative gene expression values in CTC-enriched blood fraction. Samples containing 
missing values were excluded from the analysis. On the left-hand side cluster of CTC-enrlched blood samples with elevated VIM (P = 0.068), HER2, 
CXCR4 and uPAR (P<0.00001, for all three) expression In which lymph node involvement was more frequently observed (P = 0.04). 
dol:10.1371/journal.pone.0093901.g003 



supported by research aiming at deciphering molecular mecha- 
nisms behind them. 

Conclusions 

In summary, our results show increased CTCs detection rate as 
well as more frequent VIM, SNAIL, and «Alff-positive rates of 
CTC-enriched blood fraction in patients with lymph nodes 
metastases. This indicates differences in both seeding potential 
and increased ability of the tumor cells of N+ patients to reach and 
divide in lymph nodes. We have also shown that E?M-positive 
{CK19-/VIM+) CTC-enriched blood fractions, unlike CK19+/ 
VIM— fractions, show elevated CXCR4 and uPAR levels, which 
might contribute to more aggressive tumor characteristics. 

Supporting information 

Figure SI Exemplary photos of immunostained cells 
isolated from CTC-enriched blood fractions of breast 
cancer patients. 

(PDF) 



Figure S2 Relative genes expression levels in breast 
cancer patients and controls. Expression of VIM, HER2, 
CXCR4, uPAR, TWIST! , SNAIL, CK19 and MGBl in CTCs- 
enriched blood fractions of breast cancer patients with invasive 
carcinoma (BC), patients with ductal carcinoma in situ (DCIS) and 
healthy controls (HC). As no sample expressed SLUG, results for 
this gene are not presented. 
(PDF) 

Table SI Spearman's rank correlations coefficients of 
relative gene expression levels in CTC-enriched blood 
fractions. 

(PDF) 
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